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The spectrum of gravitational waves in Randall-Sundrum braneworld cosmology 
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We study the generation and evolution of gravitational waves (tensor perturbations) in the context 
of Randall-Sundrum braneworld cosmology. We assume that the initial and final stages of the 
background cosmological model are given by de Sitter and Minkowski phases, respectively, and 
they are connected smoothly by a radiation-dominated phase. This setup allows us to discuss the 
quantum-mechanical generation of the perturbations and to see the final amplitude of the well- 
defined zero mode. Using the Wronskian formulation, we numerically compute the power spectrum 
of gravitational waves, and find that the effect of initial vacuum fluctuations in the Kaluza-Klein 
modes is subdominant, contributing not more than 10% of the total power spectrum. Thus it is 
confirmed that the damping due to the Kaluza-Klein mode generation and the enhancement due 
to the modification of the background Friedmann equation are the two dominant effects, but they 
cancel each other, leading to the same spectral tilt as the standard four-dimensional result. Kaluza- 
Klein gravitons that escape from the brane contribute to the energy density of the dark radiation 
at late times. We show that a tiny amount of the dark radiation is generated due to this process. 

PACS numbers: 04.50.+h, ll.10.Kk, 98.80.Cq 



I. INTRODUCTION 

Cosmological inflation predicts the gravitational wave 
background arising due to quantum fluctuations in the 
graviton field. Gravitational wave fluctuations are 
stretched beyond the horizon radius by rapid expansion 
during inflation, and at a later stage they come back 
inside the horizon possibly with rich information on the 
early universe and hence on high energy physics. Though 
yet undetected, gravitational waves will provide us with 
a powerful tool to probe fundamental physics in near fu- 
ture 0. 

Motivated by string theory, recently braneworld sce- 
narios have attached much attention, in which our four- 
dimensional universe is realized as a brane embedded in 
a higher dimensional bulk spacetime 0. Among them, 
the Randall-Sundrum model 0,3 is of particular inter- 
est because it includes nontrivial gravitational dynamics 
despite rather a simple construction. In the Randall- 
Sundrum type II model Q with a single brane embedded 
in an anti-de Sitter (AdS) bulk, although the fifth dimen- 
sion extends infinitely, the warped structure of the bulk 
geometry results in the recovery of four-dimensional gen- 
eral relativity on the brane at scales larger than the bulk 
curvature scale i or at low energies |fj, |6J . In order to re- 
veal five-dimensional effects particular to the braneworld 
scenario, we have to focus on the scales smaller than £, 
and for this purpose cosmol ogical p erturbations from in- 
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will be quite useful for the reason mentioned above. 

Gravitational waves from inflation on the brane were 
first studied by Langlois et al. under an assumption 
that inflation is exactly described by de Sitter space- 
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time. In this special case, the perturbation equation is 
separable and analytically solvable. A toy model called 
the "junction model" 0, is an extended version of 
the pure de Sitter braneworld, which allows a sudden 
change of the Hubble parameter H by joining two maxi- 
mally symmetric (i.e., de Sitter or Minkowski) branes at 
some time. Later, the junction model is extended to a 
more general inflation model with a smooth expansion 
rate To make the cosmological model more realis- 

tic, one should take into account the radiation-dominated 
phase that follows after inflation, and, at least in the 
low energy regime (£H <C 1), corrections to the evolu- 
tion of gravitational waves are shown to be small mm 
(see also Refs. [13, EE EE1I3)- 

In a much more general 
and interesting case, i.e., in the high energy (£H ^> 1) 
radiation-dominated phase, the perturbation equation no 
longer has a separable form and hence one cannot even 
define a "zero mode" and "Kaluza-Klein modes" without 
ambiguity. To understand the evolution of gravitational 
waves in that regime, numerical studies have been done 
by Hiramatsu et al. 0, and by Ichiki and Naka- 
mura IHEl. Their results give us a lot of implications, 
for example, on the damping nature of the gravitational 
wave amplitude due to the Kaluza-Klein mode genera- 
tion, but the initial condition they adopt is naive, ne- 
glecting initial quantum fluctuations in the Kaluza-Klein 
modes. Hence, its validity is open to question. 

The goal of the present paper is to clarify the late 
time power spectrum of gravitational waves in the 
Randall-Sundrum brane cosmology, evolving through the 
radiation-dominated stage after their generation during 
inflation. We closely follow the same line in our previous 
work |ll| , in which, using the Wronskian formulation, we 
have formulated a numerical scheme for the braneworld 
cosmological perturbations. The initial condition in our 
analysis is imposed quantum-mcchanically, and therefore 
we will be able to obtain a true picture of the genera- 
tion and evolution of gravitational perturbations in the 
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braneworld. 

This paper is organized as follows. In the next sec- 
tion, we start with giving the background cosmological 
model and summarize basic known results concerning the 
gravitational wave mode functions in the de Sitter and 
Minkowski braneworlds. In Section III, we describe the 
Wronskian formulation to obtain the power spectrum of 
gravitational waves, and then we show our numerical re- 
sults in Sec. IV. In Section V we discuss an amount of 
the dark radiation generated due to excitation of Kaluza- 
Klein modes. Finally we conclude in Sec. VI. 



II. PRELIMINARIES 

A. The background model 

Now we describe a model for the background. We shall 
work in the cosmological setting of the Randall-Sundrum 
braneworld, and so the bulk is given by a five-dimensional 
AdS spacetime. The AdS metric in the Poincore coordi- 
nates is 



ds 2 



{-dt 2 + 5 lj dx i dx j + dz 2 ), 



(1) 




de Sitter 



where £ is the bulk curvature scale and constrained by 
table-top experiments as £ < 0.1 mm |22|. A cosmological 
brane moves in this static bulk, the trajectory of which 
is given by z = z(t). The scale factor of the universe is 
related to the position of the brane as a{t) = £/z(t). 

We consider the following cosmological model on the 
brane. The initial stage of the model is given by de Sit- 
ter inflation with a constant Hubble parameter H = Hi, 
which is smoothly connected to the radiation-dominated 
phase. (In order to join the two phases smoothly, the 
brane is not exactly de Sitter at the very last stage of 
inflation.) In the radiation stage the scale factor evolves 
subject to the modified Friedmann equation |23| 



H 2 = 



1 + El 

3Af2j V 2a 



(2) 



where p r is the radiation energy density and a = 6Mpj/ ' £ 2 
is the tension of the brane. Since the conventional conser- 
vation law holds on the brane, we have p r oc a -4 . Thus, 
in terms of the proper time r on the brane we obtain 



a(r) = a(n) 



2c 



n 



2c 



1/4 



(3) 



where T\ is a fiducial time, c := ^/l + [p r (Ti)/2a] — 1, and 
p r (Ti)/a — £ 2 /8t 2 . After a period of time the energy 
scale of the universe becomes sufficiently low, and the 
radiation-dominated phase is then smoothly connected 
to the Minkowski phase. This artificial connection will 
not cause any unexpected problems on our final result 
(i.e., power spectra of gravitational waves) because at the 
end of the radiation stage the brane universe is already 



FIG. 1: Brane trajectory in static coordinates. 



in the low energy regime. Just for simplicity we assume 
that the Minkowski brane is located at z = I. Namely, 
the scale factor is normalized so that a = a = 1 when 
the universe ceases expanding. The motion of the brane 
is shown in Fig. ^ 



B. Minkowski and de Sitter braneworlds 

Let us consider tensor perturbations in AdS spacetime 
bounded by a brane. We can decompose the graviton 
field into a zero mode and Kaluza-Klein (KK) modes 
without ambiguity when the brane is maximally sym- 
metric. This is the reason why the initial and final stages 
of the background model are given by the de Sitter and 
Minkowski phases, respectively. 

We write the perturbed metric as 



ds 2 = — [-dt 2 + (5 i:j + h l3 )dx l dxi + dz 2 ] , (4) 



where hij is the transverse-traceless metric perturbation. 
As usual we decompose it into the spatial Fourier modes 
as 



h ^ 

ij (27rM 5 )V2 



d 3 k k e ik x ei 



(5) 



Here M5 is the fundamental mass scale which is related 
to the four-dimensional Planck mass Mp\ by £(M$) 3 = 



3 



Af|j. The prefactor y/2/(M 5 ) 3 / 2 is chosen so that the 
kinetic term for 0k in the effective action is canonically 
normalized. From now on we suppress the subscript k. 

For the analysis of perturbations from the Minkowski 
brane, the above Poincare coordinate system will be best 
suited, and the perturbation equation is 



dt 2 



d 2 



3 d 



dz 2 z dz 



= 0, 



subject to the boundary condition 



(6) 



(7) 



Now it is easy to find mode solutions of Eq. © . Going 
to quantum theory, the graviton field can be expanded 
in terms of the zero mode and KK modes as 

= A Q tp Q + ij^g + J dm (A m cp m + A^ipl^j , (8) 

where A n and A\ (n — 0, m) are the annihilation and 
creation operators, respectively, of their corresponding 
modes. The normalized zero mode function is given by 



<Po(t) = 



1 



-ikt 



'2ki 



while the normalized KK mode function is 

1 



<p m {t,z) = 



u m {z), 



with 

u m (z) 

and 



2 rmYi(m£)J2(rnz) — J\{m,£)Y-2,(mz) 



V[nH)] J + [JiM] ! 



= \/ k 2 + m 2 . 



(9) 



(10) 



(11) 



(12) 



The normalization here is determined by the Wronskian 
conditions 



{ifO ■ <po) = ~(Vo ■ <Po) = 1. 

(<Pm ■ <Pm') = -(<Pm ' <Pm') =S(m-m'), 

(<P0 ■ <Pm) = (<A) ' <Pm) = 0, 

(<Pn ■ ( Pn') = °) f° r n > n ' = °J m > 

where the Wronskian is defined by 

(X -Y) :=-2i dz {~^) (Xd t Y* -Y*d t X) 



(13) 



(14) 



In the de Sitter braneworld we introduce another set 
of coordinates (??,£), which is related to (t, z) as 



t = r/ cosh £ + to, z = — rysinh^, 



(15) 



where to is an arbitrary constant. In (rj, £) frame the 
de Sitter brane is located at a fixed coordinate position 



£ = & = constant, and the Hubble parameter on the 
brane is given by Hi = ^ sinh£i,. The perturbation 
equation again has a separable form 



d 2 2d 



sinh 3 £ d 1 d 



= 0,(16) 



drf -qdrj if <9£ sinh 3 £ <9£ , 

subject to the boundary condition 

M f=& = °- ( 1? ) 

Treating as an operator, the graviton held can be ex- 
panded as 

/>oo 

= a o 0o + »o0o + I dv {a v 4> u + alfil) , (18) 



where a n and ajj (n — 0, v) are the annihilation and 
creation operators of each mode. The explicit form of 
the normalized zero mode is 



Mv) = C(tHi) 



— ikrj 



(19) 



with 
C{x) 



yl + x 2 + x 2 In 



i + vi + x 2 



-1/2 



, (20) 



and the KK mode functions are found in the form of 

<f>u(r),0 = i>u(v)xAQi wh ere 

Ml) = ^- 3/2 e-^ 2 (-v) 3/2 Hll\-kr,), (21) 



P 



l/2+iv 



(cosh£) 



with 



Ci 



-C a O_? /a+4 „(coshf) 
T{iv) 



(22) 



T(h/2 + iv) 

r(-iu) 



T{h/2-iv) 



r(ii/-3/2) 



r(i+ti/) 



-i -1/2 



C*2 



P 



-1/2+t 



,(cosh£ & 



'-1/2+ 



;„(cOsh£ b ) 



,(23) 
(24) 



Note that the index v{> 0) is related to the Kaluza-Klein 
mass as 



m = I v* + - ) JZ? 



(25) 



The normalization of the modes is determined by the 
Wronskian conditions 

(00 • <t>o) = -(<t>*o ■ <t>t) = 1, (26) 

(0w • 0v) = -{4>t ■ 4>t>) = s i v - 
(^o ■ 4> v ) = (05 • 4>t) = °> 

(0„ • 0*,) = 0, for n, n' = 0, za 
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Here the Wronskian is written in (77, £) frame as 



(X ■ Y) := -2i / d£_ 



£ b rj z Sinn 4 



(27) 



III. WRONSKIAN FORMULATION 



Due to the presence of an infinite tower of Kaluza-Klcin 
modes, cosmological perturbations in the braneworld 
have infinite degrees of freedom. Instead of solving an 
initial value problem for such a system, it would be bet- 
ter to use the Wronskian formulation in order to take 
necessary degrees of freedom out of infinite information. 
In the present case, we would like to know the final ampli- 
tude of the zero mode, and therefore in fact what we need 
to do is solving the (backward) evolution of a single de- 
gree of freedom Hot lllj . Following the same line as our 
previous work |ll| . we now explain how we compute the 
amplitude of gravitational waves in the final Minkowski 
phase using the Wronskian. 

Using double null coordinates 



u = t — z, 
V = t + z, 



(28) 
(29) 



which will be convenient for numerical calculations, the 
metric Q can be rewritten in the form of 



ds 1 



U 2 



(v — u) 2 



-dudv + 5ijdx l dx- 1 ^ 



(30) 



The trajectory of the brane can be specified arbitrarily 
by 

v = q(u). (31) 
By a further coordinate transformation 



U = u, 
q(V) = v, 



(32) 
(33) 



we obtain 



ds 2 



U 2 



[q(V) - U] 



[-q'(V)dUdV + 6 i] dx 1 dx j ] , (34) 



where a prime denotes differentiation with respect to the 
argument. Now in the new coordinates the position of 
the brane is simply given by 



U = V. 



(35) 



We will use this coordinate system for actual numerical 
calculations. 

The induced metric on the brane is 



U 2 

da 2 - 

" ~ [q(V) - V} 2 



[-q'(V)dV 2 + S.jdx'dx 3 ] , (36) 



from which we can read off the scale factor a and the 
proper time r, respectively, as 



U 



q(V)-V 



dr = ayYOO dV, 



(37) 
(38) 



and hence the Hubble parameter on the brane is written 



as 



£H = 



1 



2vW) 



or equivalently 

q'{V) = (y/l+PH 2 - III}' 



(39) 



(40) 



Given the Hubble parameter as a function of r, one can 
integrate Eqs. I|38|l and (|40|l with the aid of da/dr = aH 
to obtain q as a function of V. If the Hubble parameter 
on the brane is constant in time, we have q' — constant. 
Especially, q' = 1 in the Minkowski phase. 

The Klein-Gordon-type equation for a gravitational 
wave perturbation <fi in the (U, V) coordinates reduces 
to 



(i 



uvv 



:{dv-<t(y)du)+<t{V) k 2 



q{V)-W 

supplemented by the boundary condition 



\9tj - d\ 



\u=v 



= 0. 



= 0, 
(41) 

(42) 



The expression for the Wronskian evaluated on a con- 
stant V hypersurface is given by 



(X-Y) = 2i 





2£ 


f dU 

1 — OO 


[q(V)-U\ 



(XduY*-Y*duX), 



(43) 



which is independent of the choice of the hypersurface. 

As explained in Section II, our cosmological model is 
composed of the de Sitter inflationary phase followed 
by the radiation-dominated epoch, which is connected 
smoothly to the final Minkowski phase. In the initial 
de Sitter phase, the graviton field can be expanded as 
Eq. (|18|) , while in the final Minkowski phase it can be ex- 
panded as Eq. © . We assume that initially the gravitons 
are in the de Sitter invariant vacuum state annihilated by 
So and a„, 



a |Q) = a„|0) = 0. 



(44) 



The expectation value of the squared amplitude of the 
zero mode in the final stage is 

(0|(Wo + ^4) 2 |0) = |^| 2 (0|(l + 2iti Q ) | ) 

+ oscillating part 
1 



k£ 



-„N f 



5 



where Nf :— (0|^4q Ao |0) is the number of created zero 
mode gravitons. Here we used the commutation relation 



A), A 



= 1 and assumed that Nf ^ 1. The final power 



spectrum is then given by 

ink 3 

V(k) 



(2ir) 3 (M 5 )3 
k 2 



kl 



-,N t 



(45) 



The operator Aq can be projected out by making use of 
the Wronskian relations. Noting that the Wronskian is 
constant in time, we have 

A = (<f)-(p ) f = (0-$) 

= (00 • $)ifl + J dv{4> v ■ §)ia v + h.c, (46) 

where $ is a solution of the Klein-Gordon equation l|41|) 
whose final configuration is the zero mode function ipo 
in the Minkowski phase, and subscript / and i denote 
the quantities evaluated on the final and initial hyper- 
surfaces, respectively. Thus we clearly see that final zero 
mode gravitons are created from the vacuum fluctuations 
both in the initial zero mode and in the KK modes: 



N f =\(<p*o-m 1 



dv\(4>* v ■ <$>)i 



(47) 



Correspondingly, the power spectrum (|45ll can be written 
as a sum of the two contributions: 



V = Vn 



where 



V ■= 
Vkk ~ 



tt 2 M 2 , 



TT 2 M 2 , 



Vkk, 



1(05- * 



h\ 2 , 



H«-*)i| 5 



(48) 

(49) 
(50) 



IV. SPECTRUM OF GRAVITATIONAL WAVES 

De Sitter inflation on the brane predicts the flat pri- 
mordial spectrum |8( 



4 



2C 2 (£H t 
M 2 



Hi 
2n 



(51) 



During inflation the gravitational wave perturbations are 
stretched to super-horizon scales, and then they stays 
constant until horizon reentry, with their amplitude given 
by Eq. I|51l) . (<fi = constant is a growing solution of 
Eq. (|41|l in the limit k 2 — » 0.) The primordial spec- 
trum of gravitational waves from non-de Sitter inflation 
was studied extensively in [llj, and so in this paper we 
concentrate on the simple case where inflation is given 
by the exact de Sitter model. 



For long wavelength modes with k <C fc* , where 



(52) 



labels the mode that reenters the horizon when £H = 1, 
the amplitude will decay as /i M „ cx a -1 after horizon reen- 
try, because gravity on the brane is basically described 
by four-dimensional general relativity in the low energy 
regime (£H <C 1). In fact, it is explicitly shown that 
leading order corrections to the cosmological evolution 
of gravitational waves are suppressed by £ 2 and i 2 In I at 
low energies ^201- I n particular, modes with k < fco, 
where 



(53) 



and Hq is the Hubble parameter evaluated at the end 
of the radiation-dominated phase (i.e., just before the 
Minkowski phase) 1 , reenter the "horizon" in the final 
Minkowski phase, and then they begin to oscillate (but 
their amplitude will not decay). As a result, the mean- 
square vacuum fluctuations of such modes become half 
of the initial value, leading to 



^ St 

V = — 
' 2 ' 



for k < kg. 



(54) 



In Ref. H|, the same spectrum was obtained for scales 
larger than the AdS and horizon scales by using the 
"junction model" , in which an instantaneous transition 
from a de Sitter to a Minkowski brane was assumed. For 
the reason mentioned above, modes with k < k -C fc* 
have the standard spectrum, 



(55) 



k 



Something nontrivial may happen to gravitational 
waves with k > fc*. If the effects of mode mixing are 
neglected, only the modification of the background ex- 
pansion rate alters the spectrum for these short wave- 
length modes to 



2 V k-k 



-2/3 



(56) 



However, this evaluation will not be correct because 
mode mixing is expected to be efficient at high energies. 

Our procedure to obtain a correct power spectrum is as 
follows. We solve the perturbation equation (|41|l with its 
boundary condition set to be $(E/, Vf) — (po on the final 



1 To be more precise, the scale factor in Eq. 1531 should be re- 
placed by a bit smaller one than ao = 1, because ao is defined 
as the scale factor in the Minkowski phase (where the Hubble 
parameter vanishes). However, as the smooth period connecting 
the radiation and Minkowski phases is taken to be very short, 
this point is not important. 
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hypersurface. The numerical backward evolution scheme 
we use here is the same as that used in our previous 
work yjj , and the detailed description of the scheme is 
found there. After obtaining the configuration on the 
initial hypersurface, we evaluate the Wronskian to get 
the power spectrum. 

We performed numerical calculations for three different 
values of the inflationary Hubble parameter, £Hi = 10, 
42, and 100. The radiation-dominated phase is termi- 
nated when the Hubble parameter decreases down to 
H ~ 0.03/£ =: Ho, and then it is connected smoothly to 
the Minkowski phase, so that we can see the amplitude 
of the well-defined zero mode. The numbers of grids are 
90000 in the [/-direction and 12000 in the F-direction, 
and the grid separation is chosen to be ~ 0.005 x t. In- 
tegration over the KK index v is performed up to v = 25 
with an equal grid spacing of 0.05. 

The power spectra of gravitational waves are shown in 
Fig. [3 Assuming that the power spectrum is of the form 



(a) 



V = A 



(57) 



we find that irrespective of the inflationary energy scale, 
the parameters are approximately given by 



.4 



4 



-2. 



(58) 
(59) 



Namely, we have the same spectrum as the standard one 
[Eq. I|55|) ] even for short wavelength modes with k > k*. 
As is shown in Fig. the contribution of the vacuum 
fluctuations in the initial KK modes to the final spec- 
trum, 



rKK(k) := 



Vkk 



Vq + Vkk ' 



(60) 



never exceeds 10% so far as the present calculations 
are concerned, and hence it gives a subdominant effect. 
On the other hand, the excitation of KK modes sup- 
presses the amplitude of the gravitational waves rela- 
tive to Eq. JSHJl, and our result implies that the effect 
of the modification of the background Friedmann equa- 
tion compensate this suppression, leading to approxi- 
mately the same spectral tilt as that in conventional four- 
dimensional cosmology. This is consistent with the nu- 
merical study by Hiramatsu et al. [19J, in which they 
assume the initial configuration in the bulk to be a de 
Sitter zero mode and obtain Vo oc fc~ 2 . 

Unfortunately, due to the limited number of grids in 
the [/-direction, it is difficult to evaluate accurate values 
of tkk', the convergence is not so good (Fig.^J. However, 
since tkk decreases with an increasing number of grids, it 
is strongly expected that the effect of the initial KK fluc- 
tuations is negligibly small. To obtain a more accurate 
evaluation of the contribution of the initial KK modes, we 
need an improved numerical formulation, though it seems 
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FIG. 2: Power spectra of gravitational waves from inflation 
with £Hi = 10 (a), 42 (b), and 100 (c), normalized by 8^/2. 
The total power spectrum is shown by red circles, while blue 
crosses represent the contribution only from the initial zero 
mode. Green squires indicate results from basically four- 
dimensional calculations, only including the effect of the mod- 
ification of the background expansion rate. 



quite unlikely that tkk turns to increase at a much larger 
number of grids. (We confirmed that the convergence of 
the zero mode part Vo is sufficiently good.) 



V. GENERATION OF DARK RADIATION 

So far we have concentrated on final zero mode gravi- 
tons created from initial vacuum fluctuations. In this 
section we shall discuss the generation of KK gravitons. 
In particular, we are interested in the KK mode gravi- 
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0.0001 
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FIG. 3: Contribution of initial Kaluza-Klein fluctuations, 
tkk, for the model with IHi = 100. 
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FIG. 4: r K K for IHi = 100 and k = 13.3 x fc„ versus the 
number of grids in the [/-direction. 



tons created from initial fluctuations in the zero mode, 
because from the five-dimensional point of view they are 
interpreted as gravitons that escape from the brane. At 
sufficiently late times, all the emitted gravitons fall deep 
into the bulk, and then the bulk spacetime is described 
as an AdS-Schwarzschild black hole, the mass of which 
divided by a 4 is viewed as the "dark radiation" from a 
brane observer El El El E3 • 



Before going to the estimation of the energy density 
of final KK mode gravitons, first let us take a look at 
the energy density of zero mode gravitons pew- Since 
at low energies gravitational waves evolve in a standard 
manner, their energy density behaves as Pgw <x a~ 4 - 
Therefore the ratio pgw / Pr is an invariant quantity in the 
low energy regime, irrespective of the cosmic expansion. 
Evaluating it at the end of the radiation stage, we have 



Pgw Pgw,o 



5J 
6 



(61) 



Pr Pr,0 

Note that St < 10~ 5 . In deriving the estimate (|61|l . we 



p GWa = / d\nk k 2 V(k), 



(62) 



and substituted the numerical result obtained in the pre- 
vious section, V ~ (S T /2)(k/k )~ 2 , where k can be elim- 
inated in favor of p r $ by using the Friedmann equation 
at low energies k 2 = H 2 ~ p r! o/(3Afpj). The upper limit 
of the integral may be given by the inverse horizon scale 
at the end of inflation, 



k% . CL{ H{ , 



(63) 



because the particle production is exponentially sup- 
pressed on sub-horizon scales. 

Let Pgw be the energy density of gravitational waves 
obtained by neglecting the mode mixing effect. More pre- 
cisely, pew is the energy density of gravitational waves 
h^ lL , where /i M „ is a solution of the conventional perturba- 
tion equation (d 2 +3Hd T +k 2 /a 2 )h fJil , — with the cosmic 
expansion given by a solution of the modified Friedmann 
equation. This would be much greater than pgw- Then, 
Ap := pew — Pgw — Pgw is the energy density that leaks 
from the brane, and by definition Ap is proportional to 
a~ 4 as long as it is evaluated in the low energy regime. 
Thus, Ap/p r is an invariant quantity. Now Ap ; o can be 
calculated from the spectrum of the form l|56|l . and we 
have an estimate 



Ap 

Pr 



Ap,Q 
Pr,0 
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X ffli 



(64) 



a" p r = a 4 cr = pr.o and the Friedmann 
equation at high energies, H 2 ~ p 2 /6Mp l cr. The esti- 
mate ()64|l implies that a large amount of energy (com- 
pared to pgw) is lost from the brane. Is the escaped en- 
ergy Ap directly transferred to the final bulk gravitons? 
To discuss this point, we compare it with the energy den- 
sity of the generated dark radiation. 

Since KK modes are excited dominantly at high en- 
ergies but not at low energies, the dark radiation, as is 
deduced from its name, behaves like a radiation compo- 
nent, pdr oc a~ 4 , at late times. Hence, we shall see the 
ratio Pdr/Pj- but it may be evaluated at the end of the 
radiation stage. 

The total number of the created bulk gravitons is given 
again by the Wronskian as 



dm(0\AlA m \0) = 



dm 



\{<t>*0-$m)i\ 



dv\((j>l ■ $ m ) l 



where $ m is a solution of the Klein-Gordon equation 141|) 
whose final configuration is a KK mode function (p m in 
the Minkowski phase. Concentrating on the first part, 
\{<t>t)'^m)\ 2 dm is identified as the number of KK gravitons 
coming from the initial zero mode fluctuations, with their 
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FIG. 5: Energy of bulk gravitons multiplied by the phase 
space factor. 



mass between m and m + dm. Thus the energy density 
is expressed as 



Pdr,o 



d 3 k 
(2tt) 3 

d\nk 



dm Lo\((j>l ■ $ m )j 



mk 3 uj 
d\nm „ | (0q • $ r . 



2ir 2 



(65) 



As a side remark, in the junction model of Ref. 0, the 
second part, i.e., the number of KK gravitons created 
from the initial fluctuations in KK modes, is suppressed 
compared to the first one. 

We numerically calculated the Wronskian |(0q • <5 m )| 2 
in a similar manner to the previous section. The Hub- 
ble parameter during inflation is chosen to be £Hi = 42. 
In this case, the number of grids are 150000 in the U- 
direction and 20000 in the ^-direction, and the grid sep- 
aration is about 0.0036 x £. 

The integrand 



/ 



mk 3 uj . l2 



is plotted in Fig. [SJ For each k, integration over m is 
performed up to m ~ 2/£ with a grid spacing of A In m ~ 
0.17. Performing the integration over k, we obtain 



This result indicates that the energy density of the gen- 
erated dark radiation is not larger than that of zero mode 
gravitons [Eq. (ffTTJl ]. Of course, this is a completely harm- 
less amount of an extra radiation component 28] . The 
scattering of particles on the brane in the early universe, 
discussed in [2J, [25|, |2(| |27| , can be a more efficient way 
to produce bulk gravitons. 

Although a large amount of energy is lost from the 
brane via excitation of KK modes in the high energy 
regime, the final energy density of the dark radiation is 
much smaller than that, without an enhancement fac- 
tor like £Hj in Eq. I|64|) . This discrepancy is explained 
as follows [24L l26|. In the high energy regime, the mo- 
tion of the brane is so relativistic (in the frame defined 
by the static bulk coordinates) that emitted gravitons 
run almost parallel to the brane trajectory. These gravi- 
tons stay in the vicinity of the brane and bounce off it 
many times during the high energy stage, until eventually 
they are reflected by the non-relativistic brane to fall off 
into the bulk. During this process, the gravitons lose a 
large portion of their momentum transverse to the brane 
because they repeatedly hit the retreating brane. This 
qualitatively accounts for the smallness of the final energy 
density of the dark radiation. To justify the above inter- 
pretation quantitatively, a more rigorous analysis will be 
needed in the direction of Refs. [Ig, |2^| , which includes 
calculating the pressure to the brane due to the effec- 
tive energy-momentum tensor of the bulk gravitational 
waves. 

Here we should comment on the result of the junction 
model obtained by Gorbunov et al. 0] . In terms of the 
power spectrum, their result is summarized as 



5 2 T 



2 ' 
4 



2 {k/Kf 



(k < fc*), 

(fc* < k < hi), 



(69) 



where fc* = a*/£, fc, = a^i^, and a* = a.;(= 1) be- 
cause a de Sitter inflationary stage is directly joined to 
a Minkowski phase in the junction model (see also Ap- 
pendix A of Ref. H3). From this we can estimate the 
energy density that leaks from the brane as 



p D R,o « 0.04 x r 



(66) 



where one should note that contributions from modes 
with k < k* are suppressed. The radiation energy density 
can be written as 



p r<0 = ~ a^&M^a) 1 ' 2 ^ ~ — — ( ai H t )\ (67) 



which can be obtained by using the modified Friedmann 
equation and noting that 5% ~ 3£Hf / 'in 2 Mpj for £Hi > 
1. From Eqs. I|66|l and (|67|l we have an estimate 



PDR 

Pr 



< 0(1) X 5 2 T . 



(68) 



Ap: 



M^HfS 2 T 



(70) 



where p e .i. is the energy density at the "end of inflation'' . 
On the other hand, according to Appendix D of Ref. , 
the energy density of created KK gravitons is given by 



PDR 



H 



r2 

Pe.i. Of. 



(71) 



Thus, we find that Ap ~ /?dr x £Hi, which is consis- 
tent with our present result. [Note that in the junction 
model the energy density of final zero mode gravitons 
is estimated as pew ~ (^pi/^ 2 )^t an d hence pew "C 

Pe.i.&T ~ PDR-] 
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VI. SUMMARY 

We have examined the power spectrum of the gravita- 
tional wave background in the cosmological scenario of 
the Randall-Sundrum braneworld. There are three pos- 
sible ingredients which may lead the power spectrum to 
a non-standard one: the unconventional background ex- 
pansion rate due to the p 2 term in the Friedmann equa- 
tion, the excitation of KK modes during the radiation- 
dominated stage at high energies, and the effect of initial 
vacuum fluctuations in KK modes. Previous estimates 
are based on a rather simple toy model || or numeri- 
cal studies about the classical evolution of perturbations, 
neglecting the initial KK fluctuations 0, 0|. In the 
present analysis, initial conditions are set in a quantum- 
mechanical manner and hence the effect of the initial KK 
fluctuations is included. Along the same line in Ref. [Tlj . 
we make use of the Wronskian formulation to obtain the 
final amplitude of the zero mode gravitational waves nu- 
merically. We have found that the effect of initial KK 
vacuum fluctuations are subdominant: tkk < 0.1. Our 
result confirms that the damping of the amplitude due to 
the KK mode excitation and the enhancement due to the 
modification of the background expansion rate mainly 
work, but almost cancel each other. Consequently, the 
power spectrum is basically the same as the standard 
one obtained in conventional four-dimensional cosmology. 
We believe that the cancellation between the two effects 
is a phenomenon peculiar to the radiation-dominated 



phase. To make the particularity of the radiation stage 
clear, it would be interesting to investigate consequences 
of a different equation of state parameter w(= p/p) after 
the inflationary stage. This is the next issue we plan to 
report in a future publication. 

Because of the limitation of our numerical computa- 
tion, we have not been able to give a detailed evaluation 
of how much the initial KK fluctuations contribute to the 
final power spectrum. Although it is strongly indicated 
that the initial KK effect never becomes larger than the 
present evaluation, to show that it is indeed true, we 
need to improve the numerical formulation. We will also 
report on this issue in a forthcoming paper. 

We have also estimated the energy density of the gen- 
erated dark radiation numerically, and shown that only 
a tiny amount is generated. It is smaller than the energy 
density of zero mode gravitational waves. 
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